Quantum noise thermometry for bosonic Josephson junctions in the mean field regime 
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Bosonic Josephson junctions can be realized by confining ultracold gases of bosons in multi- 
well traps, and studied theoretically with the M-site Bose-Hubbard model. We show that canonical 
equilibrium states of the M-site Bose-Hubbard model may be approximated by mixtures of coherent 
states, provided the number of atoms is large and the total energy is comparable to fcsT. Using 
this approximation, we study thermal fluctuations in bosonic Josephson junctions in the mean field 
regime. Statistical estimates of the fluctuations of relative phase and number, obtained by averaging 
over many replicates of an experiment, can be used to estimate the temperature and the tunneling 
parameter, or to test whether the experimental procedure is effectively sampling from a canonical 
thermal equilibrium ensemble. 
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I. INTRODUCTION 

Quantum degenerate Bose gases in double-well poten- 
tials exhibit coherent macroscopic tunneling dynamics, 
analogous to those in superconducting Josephson junc- 
tions [J 0, 0| • The observer can detect individual well 
populations by direct optical absorption imaging, and 
can infer the relative phase of the wave packets from in- 
terference experiments 0, [H, H|. Furthermore, atomic 
interactions can be tuned over a wide range by adjusting 
particle number and double- well parameters or by means 
of Feshbach resonances (cf. 0] for a review). 

Double-well systems are often modeled within a two- 
mode approximation by the Bose-Hubbard Hamiltonian. 
The parameters of the model are the number N ^> 1 of 
atoms, the interaction energy Eq for a pair of particles 
in the same potential well, and the tunneling coupling 
energy Ej (cf . d, 0, [l(| , which use the same notation) . 
One distinguishes the "Rabi", "Josephson", and "Fock" 
regimes [111 ] according to the relations 



Ec/Ej < N- 2 
N- 2 < Ec/Ej < 1 
1 < Ec/Ej 



(Rabi) 
(Josephson) 
(Fock) . 



Recent experiments making use of strong interactions 
deep in the Josephson regime have accomplished squeez- 
ing and macroscopic entanglement [H|. On the other 
hand, coherent tunneling dynamics and Bloch oscilla- 
tions have been realized in completely non-interacting 
Bose gases [T3, EH . The intermediate regime of moderate 
interactions (the Josephson - Rabi boundary regime) is 
virtually unexplored in experiment. This regime, where 
N 2 Ec ~ Ej, is of particular interest, as it contains most 
of the stationary Josephson modes, such as and tt phase 
modes, and the onset of macroscopic quantum self trap- 
ping [14| . Here, number and phase fluctuations are sen- 
sitive to the ratio of N 2 Ec to Ej, both in the ground 
state [TEj and, as we shall see, in thermal equilibrium at 
higher temperatures. 



When a gas of ultracold bosons is released from a 
double- well potential trap and recombined in free expan- 
sion, interference fringes, analagous to those of Young's 
double-slit experiment, are observed in the atomic den- 
sity. This does not necessarily mean that the double-well 
system was prepared in a coherent state, for individual 
images or "shots" will feature interference fringes even if 
the gases are initally independent [16]. To ascertain that 
the experimental procedure prepares the system in a co- 
herent state, one needs to repeat the experiment many 
times and compare the results. If the fringes always lie 
in the same position, one may infer that the state is co- 
herent, and ascribe a definite value to the relative phase 
between the condensates in the two wells. 

The decohering effect of temperature is seen in the fluc- 
tuations, from one shot to another, of the location of the 
interference fringes. These fluctuations reduce the visibil- 
ity of the interference fringes when the density profiles are 
averaged. The fringe contrast in the average density pro- 
file is called the "coherence factor" and, for double-well 
systems in the Josephson regime, is found to be a certain 
function of k B T/Ej @. This function is used to calibrate 
the "thermometer" of noise thermometry [l(| [ItJ ■ 

In this article we study canonical thermal equilibrium 
states of N 3> 1 bosons in double-well and multi-well 
potentials, focusing on regimes where both Ej/k B T <C 
N and N 2 E c /k B T < N. This includes the Rabi- 
Josephson boundary regime, provided the temperature 
is high enough that Ej/k B T < N. We find that 
the coherence factor is sensitive to the ratios Ej/k B T 
and N 2 Ec/k B T. Our results are rigorous inasmuch as 
they are derived from a general theorem about canonical 
statistics of M-mode boson models [l8j |. 

The regimes we consider are not normally attained 
in atom interferometry experiments. For example, the 
noise thermometry experiments reported in [Tol . [rn ] in- 
volved trapping a few thousand 87 Rb atoms at 15-80 
nK. The parameter Ej/k B T ranged between about 0.02 
and 20, but the parameter N 2 Ec/k B T was large be- 
cause the experiments were performed deep within the 
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Josephson regime. However, it should be possible to engi- 
neer bosonic Josephson junctions in the Rabi- Josephson 
boundary regime by taking advantage of Feshbach reso- 
nances to reduce the interaction parameter Eq [HI, Ell • 
This paper is organized as follows. In Section [IT] we 
state our main results and proposals concerning noise 
thermometry of 2-site bosonic Josephson junctions. In 
Section [TTT1 we state the central result of this article, The- 
orem[T] a general result concerning canonical statistics of 
Af-site Bose-Hubbard models. We return to the particu- 
lar case M — 2 in Section HVl and discuss the fluctuations 
of density observables such as interference fringes in time- 
of-flight matter wave interferometry. We outline a proof 
of Theorem Q] in the Appendix. 



II. NOISE THERMOMETRY WITH BOSONIC 
JOSEPHSON JUNCTIONS 

Double-well s yste ms may be modeled using a two-mode 
approximation [l9(. In a symmetric two- well potential, 
the ground state is degenerate when the wells are sepa- 
rated by an infinitely high barrier: the gerade and unger- 
ade modes have the same energy. If the barrier between 
the wells is finite, tunneling lifts the ground state de- 
generacy. Provided the tunneling barrier is not too low, 
the energy splitting of these low energy modes is small 
compared to the energy difference between them and the 
higher excited states, and, at low enough temperatures, 
there are effectively only two modes in play. 

Two-mode approximations have been derived using ei- 
ther a "semiclassical" or "second-quantized" approach. 
The former approach attempts to constrain the semi- 
classical Gross-Pitaevskii dynamics to a two-dimensional 
subspace (T], [l4| . The latter approach begins with the sec- 
ond quantized Hamiltonian, and attempts to restrict it 
to one involving only two modes ppj , which might be the 
lowest energy solutions of the Gross-Pitaevskii equation 
or which might be found by self-consistent varia- 
tional minimization of the energy over all suitable two 
mode approximations [2l| . The relationship between the 
second-quantized and the semiclassical theories is dis- 
cussed in Q. 

Two-mode models become quite sophisticated [22j]. 
The simplest one is the 2-site Bose-Hubbard model, 
whose Hamiltonian is 



H 



N 



-Ejx + (A 2 £ c /4) z 2 



(1) 



In this formula, the operators x — -^{a\a2 + a\ai) and 

z = jj{a\ai — a\a2) are understood to operate on the 
iV-particle component of the boson Fock space. The ob- 
servable z is the relative number imbalance between the 
wells. The observable x is the relative occupation differ- 
ence of the gerade and ungerade modes, which is related 
to relative phase. 

Canonical thermal equilibrium states of the 2-site 
Bose-Hubbard Hamiltonian can be approximated by cer- 
tain mixtures of coherent states. We shall show that this 



approximation is rigorously justified for regimes where 
Ej/k B T < N and N 2 E c /k B T < N (and N is large). 
Define the dimensionless parameters 

5 = Ej/k B T 

e = N 2 E c /4:k B T . 

For regimes where <5, e <C N , we will derive the follow- 
ing formulas for the coherence factor ((2]) and the second 
moment of z (|5j). 

The coherence factor a is defined to be the fringe con- 
trast in the ensemble averaged density profile of a double- 
well interference experiment 0, [I?], El- We will show 
that 



/^ 1 xJ (e(l - x 2 )/A)e 5x+ex2 / A dx 
/^ 1 / (e(l - x 2 )/A)e Sx +^ 2 / 4 dx ' 



(2) 



where Iq denotes the modified Bessel function of the first 
kind (of order zero). In the non-interacting case, when 
e = 0, formula @ reduces to 



Q!(5,o = coth(<5) — 1/(5 



(3) 



In the strongly repulsive case, when s ^> 1, the term 
Jo(e(l — x 2 )/A)e ex / 4 is nearly proportional to 1/yT — x 2 
over much of the domain of integration, and formula (|2|) 
tells us that 



as. 



V 1— 37 



--dx 



V 1— ar 



-dx 



h{5) 
Io(5) 



(4) 



This agrees with the semiclassical formula for the coher- 
ence factor in the Josephson regime @, [13, [13] • 

In a symmetric double-well potential, the expected 
value of the population imbalance z is zero, i.e., (z) = 0. 
We will show that the variance of z is 



(z 2 ) = 



^z 2 h(5VT^z 2 ~)e-" ' 2 dz 
ft I {SVT^)e-" 2 /2dz 



(5) 



In Figures 1 and 2, a and (z 2 ) are plotted against 
k B T/Ej on a logarithmic scale, for various values of 



A 



e/S = N 2 E c /4Ej 



The method of "noise thermometry" developed in 
[l(J E3] uses statistical estimates of a, obtained by repli- 
cating a double-well experiment under identical condi- 
tions, to determine k B T/Ej. If Ej is known, or es- 
timable, the temperature T can be deduced, even when 
this temperature is so low that it cannot be found by the 
usual technique (fitting a gaussian to the "wings" of the 
density profile after some time-of-flight). This method is 
suitable for the Josephson regime Ej -C N 2 Ec, where 
the parameter A is formally equal to oo. 

To perform noise thermometry in the Rabi-Josephson 
boundary regime, where 0< A<C oo, one needs to know 



3 



i 

0.9 ■ 
0.8 ■ 
0.7- 

s 

2 0.6 ■ 

g 0.5 ■ 
el 

CD 

£ OA- 

o 

0.3 ■ 
0.2- 
0.1 ■ 
- 



Rabi regime (A = 0) 
•■A= 1 
-A = 5 
-A= 10 

- Josephson regime (A = ») 



10 



10" 10' 
temperature [k B T / E ] 



FIG. 1: Coherence factor a vs. temperature for several values 



N 2 Ec/AEj (logarithmic scale). 
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FIG. 2: Variance of the normalized population imbalance (z 2 ) 
vs. temperature for various values of A = N 2 Ec/iEj. The 
gray dashed lines indicate the "low temperature" approxima- 
tion (fT5j) . 



both UbT/Ej and A. Estimates of ksT/Ej and A can 
be deduced from statistical estimates of the coherence 
factor a and the variance (z 2 ) of the number fluctuations, 
thanks to formulas ((2]) and (J5]). 

There are a couple of benefits of doing noise thermom- 
etry in the Rabi- Josephson boundary regime: 

1. When a bosonic Josephson junction is fashioned in 
the laboratory, one usually has more accurate knowledge 
of the parameter Eq than the parameter Ej (the tun- 
neling energy Ej is quite difficult to estimate, due to 
its exponential sensitivity to the precise geometry of the 
double- well potential). By performing noise thermom- 
etry in the Rabi-Josephson boundary regime, one can 
take advantage of one's knowledge of Eq to estimate Ej 
as well as T. 

2. If one does happen to know Ej with some accuracy, 



one obtains two estimates of the temperature. If these 
estimates differ significantly, it may be evidence that the 
experimental procedure has failed to prepare the double- 
well system in a canonical thermal equilibrium state. The 
assumption that replication of the experiment is sam- 
pling from the canonical ensemble ought to be tested, 
because some ways of preparing the system can fix it 
in a non-canonical equilbrium state, for example, if the 
double-well is formed by ramping up a potential barrier 
too quickly (23j. 



III. CANONICAL STATISTICS OF THE M-SITE 
BOSE-HUBBARD MODEL 

The M-site Bose-Hubbard Hamiltonian for bosons 
with nearest-neighbor hopping is 



J 



M-l 

E 

=i 



(a\a l+1 +a] +1 ai) + — 



i=l 



a]a]a,iai 



We are going to discuss systems of exactly N bosons, and 
take a limit N — > oo. Accordingly, we will consider the 
restriction of the above Hamiltonian to the A^-particle 
subspaccs of the boson Fock space, and allow the pa- 
rameters J and U to depend on N. Let Vn denote the 
orthogonal projector onto the A^-particle component of 
the boson Fock space over C M , and let 



JU-l 



— « U M 
- Jn ^ (a\a i+ i + a\ +1 ai) + V,,',, 



2 ^ 

i=l 



be the the iV-boson Hamiltonian H^. 
The density operator 



alalaia^jVN 
(6) 



e -H N /k B T /Tr(e~ Hjv/fcBT ) 

represents the canonical ensemble of TV bosons in thermal 
equilibrium at temperature T, in the sense that 

(X) NtT = TY(Xe- H "/ kBT )/Tr(e- H "/ k * T ) (7) 

is the expected value of any observable X when the sys- 
tem is in the canonical thermal equilibrium state. We 
are going to show that this state may be approximated 
by a mixture of coherent states, provided that N 3> 1 
and both NJ N /k B T, N 2 U N /k B T < N. 

We parameterize the pure states of the M-site system 
by the product of the standard (M — l)-dimensional sim- 
plex and the Af-dimensional torus. Let Am denote the 
standard (M — l)-dimensional simplex 

[p= (pi,P2,---,Pm) e R M \ Pl + ---+ PM = 1 , Pi>0} 

and let [0, 2tt) m denote 

(0i,(/> 2 ,...,<M| <k € [0,2vr) for i = 1,2,...,m}. 



4 



To each point (p, <f>) € Am x [0, 2tt) m we associate the 
unit vector 

</W = (VpIe #1 ,V^e*,...,^e^) . (8) 

The parameterization (p, <f) i— > ip Pl ^ is many-one, because 
a global change of phase in <f> does not change ip Pj< j,. 

Let fiM (dp) denote the uniform probability measure on 
Am; in particular, /i 2 is equivalent to the length measure 
dp on the unit interval p € [0, 1]. Let UM{d(f>) denote the 
uniform probability measure (2ir)~ M d<fiid4>2 • • • d<pM on 
[0,2^) M . 



Theorem 1 Let J/v and Un be two sequences of param- 
eter values. For each TV , let denote the operator 
and let (X)n,t denote the ensemble average ([?[) for the 
canonical ensemble at temperature T . 

If 

lim - — ^- = S and lira — ; — = e , (9) 

iV-oo k B T N^oo Ak B T y ' 

then, for any vectors \i, ■ ■ ■ Xk, Xi, ■ ■ ■ , Xk e cM > 



lira -^r(at,aL 



l X k /N,T 



A M J[0,2tt) 



J]_(Xt\tpp,<t>)(tpp- 



exp 



M-l 

2( 5 J2 y/ViPi+1 COS(0 1+ 1 



M 
i=l 



J J e 2& T, ^Jp'iP' i+1 cos(^ +1 -<^) 



-#)- £ I>i : 



fXM(dp')l>M(d<f>'j) 



VM(d(j))fJLM{dp) 



(10) 



r 



This theorem can be deduced from propositions con- 
cerning Finetti representations for canonical states of M- 
mode bosons [18( . A proof is outlined in the appendix. 



IV. DERIVATION OF FORMULAS © AND © 



The 2-site Bose-Hubbard Hamiltonian is [8J| 

E E 
- — {a\a 2 J ra\a 1 )'P N + — (a\a\aiai + a\a\a 2 a 2 )V N 

(Vn restricts the operators to the TV-particle component 
of the Fock space). When expressed in terms of the ob- 
servables 



= jf(a\a 2 + a\ai)V N 
= jf{a[ai - a\a 2 )VN 



on the TV-boson space, this Hamiltonian only differs by 
the constant E C N(N ~ 2)/4 from the Hamiltonian (Tj). 

We are going to rewrite formula (fT0)l for M = 2, the 
double-well case. Changing variables 



pi = \ + \z, p 2 = \ 



\z, 4> = 4> 2 



in formula ([! 
Define 



we write 



n+z -i<i>/2 

V2 6 



V2 



{T±z -i<t>/2 yT^I icj>/2\ 



V2 



V2 



(11) 



for (z,(j>) G [—1)1] x [0,2tt). The operators a\ and a\ 
in Hn are identified with the creation operators for the 
vectors ui.o = (1,0) and U—ifi = (0, 1), respectively. Let 
us also define the probability density functions 



Pi,e( z >' 



exp (<$Vl — z 2 cos 4> — ez 2 /2) 

J J e -5\/l-2' 2 co S 0'-«' 2 /2^/^ z / 



(12) 



on [—1,1] x [0, 27r). Changing variables in formula ijTDl 
we find that 



lira — r(at 

2tt r l k 



■ a L a x[ 



•(-On,! 



Y[ (Xi I «* , 4, ) (u z , ^ | x- ) P 5)E ( z , 0) dzd^ 
»=i 

(13) 

in the limit TV — > oo with 

Ej/k B T — » J, N 2 E c /±k B T — » e . (14) 



A. Population imbalance 

Formula (fT5|) may be applied directly to the observable 
z. In a symmetric double-well, (z) = 0. Higher moments 
of z are those of the probability distribution 



/ (<5vT 



Jo J (5 VI - (^') 2 )e _£(z ' )2/2 ^' 



<iz 



5 



that is, in the limit (|T4[) for each fixed k, 

f\ z k I (5VT^)e-" 2 / 2 dz 



lim /; '' 



z )n,t 



f\l (6VT^)e-^ 2 / 2 dz 
We demonstrate this for k = 2: 



lim (z 2 ) 



N.T 



lim 
lim 



1 
1 

7V2 



((ajai - 4a 2 ) 2 ) 



N.T 



I T2 2 , t2 2 ott \ 

(^aj a x + 0,2 a 2 — 2a{a2aia2y 



N.T 



lim + 4 a 2) 



2tt /.l 



N.T 



{Ul,0\Uz,4>)\ - \\U-l.o\Uz,. 

i-1 

X P,5 je (2:, </>) (fe<i</> 



z 2 Ps >E (z,(j>) dzdcj) 



f^ 1 z 2 I (SVl^ 1 )e-" 2 / 2 dz 
j 1 _ 1 I (5y/l^z 2 )e-" 2 / 2 dz 



This proves formula ([5]) for 



z 2 )s, e - Figure [2] shows that 
(15) 



5 + e 



is a good approximation at lower temperatures. 



B. Coherence factor 

In a time-of-flight (TOF) experiment on double-wells, 
the potential trap is suddenly shut off and the gas ex- 
pands into free space for awhile before it is imaged. The 
images constitute a measurement of the "integrated den- 
sity" observable J (r)a(r)dr3, where a^(r) denotes the 
usual field operator at r = (n.,^,^), and the integral 
is over the spatial coordinate parallel to the imag- 
ing light beam and perpendicular to the line that passes 
through the two wells, the ri-axis. We turn our atten- 
tion to the density observables a^(r)a(r) and their linear 
combinations. 

Moments of such density observables are easily com- 
puted if atom-atom interactions during the TOF are 
neglected. In a two-mode approximation, each vector 
u z ,4> S C 2 is identified with some wavefunction ^ Ztt p(r). 
The vectors u-i ; o and ui.q are identified with the wave- 
functions of the "left" and "right" wells, respectively. 
The specific map depends on the two-mode ap- 

proximation adopted, but the precise form of the initial 
left and right well wavefunctions hardly affects the in- 
terference pattern observed after a long TOF, and we 
may simply assume that $-1,0 and \&+i,o are gaussian 
wave packets centered at (-d/2, 0, 0) and (d/2, 0,0) (2i |. 
After a long enough [25| time t of free expansion, the 



wavefunction ^±10 (r, t), which describes the state of an 
atom that was initially in right (+) or left (-) well, will 
be nearly proportional to 



exp 



d 2 /4))exp( Ti^~ t n 



(16) 



over the region where the density is imaged. Let 
^2,0(1", t) denote 

Hle-^ 2 vl/_ li0 (r,t) + ^e^H +h0 (r,t) . 



V2 



Supposing that atom-atom interactions during the period 
of expansion may be neglected, the state of the many- 
boson system at time t is just the one freely induced by 
the 1-particle map u z .$ 1— > ^ z ,<f>> an d Theorem Q] implies 
that 

1 



r2ir r l k 



J™L Tfk ( flt ( r i) ■ • ■ a t (r fc )a(r'i) ■ • • a(r' k )) N T 

/ / Jl x S'z.4,(r l ,t)'¥ z ~^(r' i ,t)Ps, £ {z,(j))dzd<f) 
Jo J-i i=1 

(17) 

for all k and points ri, . . . , and r'i, . . . , r'^, in the limit 
(fT4| . Substituting (fTBj) into (jTTJ) and proceeding formally, 
one finds that 



lim — (aUr)a(r)) 



N,T 



2tt 



oc 

where 

OiS,s 



% z ,<j,(r,i)\ P s , £ (z,(/))dzd<p 
(md 



1 

1 /-27T 



y/l-z- 
'-1 JO 

(md \ 



V ht 



r% — (p\Ps l£ (z, 4>)d<pdz 
(18) 



2t: 



\J\ — z 2 cos 4> Ps <£ (z,<p)d(pdz . (19) 



1 Jo 



Finally, fomula |(5J) for a$ t£ is obtained by changing vari- 
ables x = Vl — z 2 cos<ft, y = yl — z 2 sin^> and integrat- 
ing over y. 

Formula (fT5|) shows that the interference pattern will 
feature fringes with spacing ht/md and contrast equal 
to the coherence factor a$ i£ . In particular, formula (fl8)) 
implies that 

lim 1/ f a t( r ) a (r)e jk - r dr\ 

JV^oo N\J I N,T 

is proportional to as, e when k = (md/ht, 0, 0). Thus 
the coherence factor can be estimated by averaging, over 
many replicates of a TOF experiment, the Fourier coef- 
ficient of the imaged density profiles at wave-vector k. 
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V. CONCLUSION 

We have studied the canonical statistics of phase and 
number in the M-site Bose-Hubbard model ©. Theo- 
rem [J provides a convenient way to approximate canon- 
ical thermal equilibrium states by mixtures of coherent 
states. From Theorem [T] we have derived formulas ^ 
and (|5|) for the coherence factor a and the variance of the 
relative population imbalance in symmetric double-well 
bosonic Josephson junctions. These formulas are valid in 
the Rabi- Josephson boundary regime, provided N 1 
and Ej/k B T < N. 

We have proposed a way to perform noise thermometry 
in the Rabi- Josephson boundary regime. In this regime, 



canonical statistics depend on two parameters, e.g., the 
dimensionless parameters Ej/ksT and A = N 2 Ec/iEj. 
Statistical estimates of the coherence factor and the vari- 
ance of the number fluctuations can be used to obtain em- 
pirical estimates of Ej and T, and to test the assumption 
that the system is being prepared in a canonical equilib- 
rium state. 



Acknowledgment: A. G. is supported by the Vienna 
Science and Technology Fund project "Correlation in 
Quantum Systems" . This work was done under the 
auspices of Joerg Schmiedmayer's Atom Chip Lab. We 
thank Igor Mazets for helpful comments. 



[1] A. S. S. Fantoni, S. Giovanazzi, and S. R. Shenoy, Phys. 

Rev. Lett. 79, 4950 (1997). 
[2] R. Gati, M. Albiez, J. Foiling, B. Hemmerling, and M. K. 

Oberthaler, Appl. Phys. B 82, 207 (2006). 
[3] S. Levy, E. Lahoud, I. Shomroni, and J. Steinhauer, Na- 
ture 449, 579 (2007). 
[4] J. Esteve, C. Gross, A. Weller, S. Giovanazzi, and 

M. Oberthaler, Nature p. doi:10.1038 (2008). 
[5] Y. Shin, M. Saba, T. Pasquini, W. Ketterle, D. Pritchard, 

and A. Leanhardt, Phys. Rev. Lett. 92, 050405 (2004). 
[6] T. Schumm, S. Hofferberth, L. M. Andersson, S. Wil- 

dermuth, S. Groth, I. Bar- Joseph, J. Schmiedmayer, and 

P. Kriiger, Nature Phys. 1, 57 (2005). 
[7] T. Kohler, K. Goral, and P. S. Julienne, Rev. Mod. Phys. 

78, 1311 (2006). 
[8] J. R. Anglin, P. Drummond, and A. Smerzi, Phys. Rev. 

A 64, 063605 (2001). 
[9] L. Pitaevskii and S. Stringari, Phys. Rev. Lett. 87, 

180402 (2001). 

[10] R. Gati, B. Hemmerling, J. Foiling, M. Albiez, and M. K. 

Oberthaler, Phys. Rev. Lett. 96, 130404 (2006). 
[11] A. J. Leggett, Rev. Mod. Phys. 73, 307 (2001). 
[12] M. Gustavsson, E. Haller, M. J. Mark, J. Danzl, 

G. Rojas-Kopeinig, and H.-C. Nagerl, Phys. Rev. Lett. 

100, 080404 (2008). 
[13] M. Fattori, C. D'Errico, G. Roati, M. Zaccanti, M. Jona- 

Lasinio, M. Modugno, M. Inguscio, and G. Modugno, 

Phys. Rev. Lett. 100, 08040 (2008). 
[14] S. Raghavan, A. Smerzi, S. Fantoni, and S. R. Shenoy, 

Phys. Rev. A 59, 620 (1999). 
[15] J. Javanainen and M. Y. Ivanov, Phys. Rev. A 60, 2351 

(1999). 

[16] M. R. Andrews, C. G. Townsend, H.-J. Miesner, D. S. 
Durfee, D. M. Kurn, and W. Ketterle, Science 275, 637 
(1997). 

[17] R. Gati, J. Esteve, B. Hemmerling, T. B. Ottenstein, 
J. Appmeier, and A. W. and. M K Oberthaler, New J. 
Phys. 8, 189 (2006). 

[18] A. D. Gottlieb, Journal of Statistical Physics 121, 497 
(2005). 



[19] J. Javanainen, Phys. Rev. Lett. 57, 3164 (1986). 

[20] G. J. Milburn, J. Corney, E. M. Wright, and D. F. Walls, 

Phys. Rev. A 55, 4318 (1997). 
[21] R. W. Spekkens and J. E. Sipe, Phys. Rev. A 59, 3868 

(1998). 

[22] D. Ananikian and T. Bergeman, Phys. Rev. A 73, 013604 

(2006) . 

[23] R. Gati and M. K. Oberthaler, J. Phys. B 40, R61 

(2007) . 

[24] A. Imambekov, V. Gritsev, and E. Demler, arXivxond- 

mat/0703766vl (2007) . 
[25] t so large that ^Jht/m is much greater than the width of 

the wells. 



VI. APPENDIX: PROOF OF THEOREM [T] 

Recall that Vn denotes the projector whose range is 
the A'-particle component of the boson Fock space over 
C , and recall the notation introduced around formula 
©. Proposition 2 of [HI implies that 

n™oo ^ Tr ^ a xi ■ ■ ■ °L a xi ' ' ' a x ' fc )^)/ Tr (^) 

= 11 TT(XilVw)WwlXi) v M{d(j))^M{dp) 

(20) 

for any vectors Xii • ■ ■ Xk, Xi, ■ • • i Xk ^ C M . Indeed, for- 
mula (PO)) holds even if the product of the a^. and a x > is 
not normally ordered. 

Writing the operator defined in ^ as 

H N = ( - J N X! + (U N /2)X 2 )V N , 
we may write 



7 



oo oo 

-H N /k B T 



c 

m —0 n 2 — 



I 

We are going to take a limit of the trace of both sides According to formula (f20|) , the limit in (|2T|) equals 
of the preceding equation. Formula (I20[) implies that a 
limit such as 



is equal to the same limit for the normally ordered form 
ie operator, i.e., the limit here is identical to 

N n! + 2nJ T ( ■■ X i 1X 2 2 : V N )/Tk{V N ) . (21) g(p, 0) = £ ?l Therefore, 



M-l 

of the operator, i.e., the limit here is identical to where f(p,<fi) = 2 cos (^i+i — 0i) an d 

i=l 

1 M 



i=l 



Tr(e-^/^) = g g H>l r£ »{^^ / / fip^Mr ,M(#)MM(rfp) 

= / / E E ^w(^)) ni (^j(p^))" ! VMm»M{d P ) 



e 8t(p,*)-eg(p,4,) VM (d<l>)li M (dp) 



Am J[0,2iry 

Similarly, 

Tr((at • • -at a Y - ■ • •a Y 0e" ffjY / fcBT ) /" /" A , ^ , ^ 

lim X " T Xl ^ = / / H(XiW P ,*W P M> Sf ^- £aM) v M {d<t>)n M {dp) 

N ^°° tr{l J N) J a m J[0M) m f=l 



The last two equations imply formula (fT0|) . 



